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1.0  INTRODUCTION 


1.1  GENERAL 

The  formation  of  ice  on  aircraft  surfaces  occurs  during  flight  through  clouds  of 
supercooled  water  droplets.  Ice  accretion  on  these  surfaces  usually  results  in  a  degradation 
of  both  aircraft  performance  and  operational  safety.  Therefore,  protective  devices  for 
removing  accumulated  ice  (deicing  devices)  or  for  continuously  maintaining  the  exposed 
surfaces  of  the  aircraft  free  of  ice  (antiicing  devices)  are  routinely  used.  For  safety  reasons, 
the  effectiveness  of  the  protective  devices  must  be  determined  before  their  use;  hence,  each 
device  must  be  tested  under  simulated  atmospheric  icing  conditions. 

Icing  conditions  are  simulated  by  duplicating  the  principal  factors  that  characterize  an 
icing  cloud:  (1)  air  temperature,  (2)  droplet  size  distribution  (including  the  “mean  effective” 
droplet  diameter  for  the  droplet  size  distribution),  and  (3)  cloud  liquid  water  content 
(LWC),  The  natural  icing  conditions  simulated  during  testing  are  well  known  and 
documented  in  the  “Airworthiness  Standards”  (Ref.  1).  Figure  1  summarizes  the 
meteorological  conditions  used  for  icing  tests. 

During  icing  cloud  simulation,  the  air  temperature  is  measured  and  controlled  by 
conventional  methods  (Ref.  2).  However,  measurement  and  control  of  the  cloud  droplet  size 
distribution  and  LWC  require  specialized  instrumentation  systems  and  test  techniques.  To 
produce  true  icing  simulation  requires  careful  consideration  of  the  test  techniques  and  of 
certain  fundamental  limitations  of  the  droplet  sizing  instrumentation  systems. 

1.2  BACKGROUND 

The  engine  test  cells  at  the  Arnold  Engineering  Development  Center  (AEDC)  Engine 
Test  Facility  (ETF)  have  been  developed  to  provide  water  droplet  fields  that  simulate  cloud 
condensate  for  engine  icing  tests.  These  facilities  were  developed  to  determine  the 
operational  effectiveness  of  both  turbojet  engines  and  their  icing  protection  systems  under 
conditions  that  closely  simulate  or  duplicate  natural  icing  conditions. 

Icing  clouds  are  produced  at  the  AEDC/ETF  by  injecting  a  continuous  spray  of  water 
droplets  into  a  cold  airstream  directed  at  the  engine.  Figure  2  shows  that  the  water  spray  is 
injected  by  a  bank  of  pneumatic  spray  nozzles  located  in  the  low-velocity  region  of  the  test 
cell  upstream  of  the  test  cell  bellmouth.  The  LWC  and  the  droplet  size  distribution  (from 
which  the  mean  effective  droplet  diameter  can  be  obtained)  are  maintained  through 
variations  in  the  water  and  air  supply  pressures  of  the  pneumatic  spray  nozzles. 
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The  icing  cloud  LWC  and  mean  effective  droplet  diameter  are  specified  before  testing  as 
part  of  the  test  matrix;  normally,  more  than  one  icing  cloud  are  simulated  during  a  given 
test.  Consequently,  the  spray  nozzles  are  calibrated  before  testing  to  yield  LWC  and  mean 
effective  droplet  size  as  functions  of  the  inlet  conditions  to  the  spray  nozzles.  Although  the 
spray  nozzle  calibration  is  used  during  testing  to  establish  the  correct  icing  conditions,  it  is 
desirable  to  determine  continuously  the  accuracy  of  the  cloud  simulation.  Therefore,  the 
LWC  and  droplet  size  distribution  are  monitored  throughout  the  actual  icing  tests. 

To  characterize  the  LWC  and  mean  effective  droplet  size,  the  AEDC/ETF  uses  the  total 
injected  water  flow  rate  and  holograms  of  the  droplet  cloud  taken  at  the  engine  inlet.  The 
LWC  is  determined  from  the  measurement  of  total  injected  water  flow  rate  by  an  icing 
simulation  computer  program  described  in  Ref.  3.  The  mean  effective  droplet  size  is 
determined  from  the  reconstructed  holograms. 

It  should  be  noted  that  use  of  the  holocamera  to  determine  droplet  size  distributions  does 
not  provide  the  desired  online  continuous  review  of  droplet  information  because  the 
holographic  data  are  reduced  off  line.  However,  it  is  reported  that  several  particle  diagnostic 
systems  have  been  developed  that  can  size  particles  continuously  and  in  nearly  real  time 
(Refs.  4  through  12). 

It  was  decided  that  because  of  the  need  to  monitor  and  set  droplet  size  during  icing 
simulation,  in-depth  evaluation  of  particle-sizing  techniques  and  systems  was  required. 

1.3  SCOPE  OF  INVESTIGATION 

A  research  program  was  conducted  at  the  AEDC/ETF  to  refine  the  measurement 
techniques  used  to  characterize  the  injected  water  droplets  of  a  simulated  icing  cloud.  The 
objectives  of  the  study  were  (1)  to  develop  improved  measuring  techniques  for  defining 
water  droplet  number,  size,  and  velocity  and,  (2)  if  possible,  to  use  these  measurement 
techniques  to  obtain  a  data  base  for  an  improved  icing  simulation  computer  program 
(Ref.  3). 

To  achieve  the  program  objectives  required  (1)  surveying  available  particle  diagnostic 
systems  and  selecting  the  systems  meeting  the  measurement  requirements  defined  in  Section 
2.2,  (2)  experimentally  evaluating  these  systems  to  determine  their  applicability  as  particle 
diagnostic  tools  for  icing  simulation,  and  (3)  developing  a  research  test  cell  to  provide  the 
required  icing  environments  for  the  experiments.  The  selected  panicle  diagnostic  systems  are 
described,  and  the  droplet  size  data  obtained  during  the  experimental  program  are 
presented.  The  research  test  cell  is  also  described  and  its  operating  characteristics 
summarized. 
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2.0  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

2.1  ICING  RESEARCH  TEST  CELL 

A  research  test  cell  has  been  developed  for  the  sole  purpose  of  studying  icing 
environments.  The  research  cell  -  a  subscale  version  of  the  AEDC  propulsion  development 
test  cells  -  can  provide  the  meteorological  icing  condition  described  by  Fig.  1  for  airflows  up 
to  30  lbm/sec.  Table  1  summarizes  the  research  test  cell  characteristics  pertinent  to  the  icing 
experiments. 

The  icing  research  test  cell  (Figs.  3a  and  b)  consists  of  a  flowmetering  venturi,  plenum 
chamber,  water  spray  system,  bellmouth,  removable  inlet  connecting  ducts,  and  a  test 
chamber.  A  secondary  air  supply  system  also  provides  air  to  the  test  chamber;  this  secondary 
air  encapsulates  the  primary  air  to  prevent  recirculation  of  the  icing  cloud  water  droplets 
through  the  measurement  plane  (Fig.  3a). 

The  plenum  chamber  rides  on  casters  along  rails.  It  can  be  moved  in  the  axial  direction; 
hence,  the  axial  distance  between  the  spray  nozzle  (located  in  the  plenum  chamber)  and  the 
test  section  can  be  varied  without  moving  the  test  chamber.  This  makes  it  possible  to  study 
water  droplet  formation  as  a  function  of  axial  position  without  relocating  the  particle 
diagnostics  instrumentation. 

The  water  spray  system,  which  provides  filtered,  demineralized  water  and  atomizing  air 
to  the  spray  nozzle,  is  shown  in  Fig.  3c.  A  single  Spraying  System  26B  spray  nozzle  (Fig.  3d) 
provides  the  droplet  cloud.  The  spray  nozzle  is  located  on  the  centerline  of  the  plenum 
chamber  and  is  attached  to  a  movable  screw  assembly.  The  axial  position  of  the  spray  nozzle 
can  be  varied  from  12  to  18  in.  upstream  of  the  bellmouth. 

The  test  chamber  (Figs.  3a  and  b)  has  six  instrumentation  and/or  viewing  ports  around 
its  circumference.  The  icing  cloud,  which  issues  from  the  connecting  duct  into  the  test 
chamber  as  a  free  jet ,  can  be  viewed  through  all  of  the  ports.  The  test  plane  -  the  plane  along 
which  all  droplet  size  measurements  are  made  -  coincides  with  the  centerline  of  the  ports. 

2.2  ICING  INSTRUMENTATION 

For  a  particle  diagnostics  system  (PDS)  to  meet  the  testing  needs  of  the  AEDC/ETF,  it 
must  be  capable  of  determining  particle  diameters  and  number  densities*  in  the  range 


*For  a  given  LWC  in  Fig.  1,  the  droplet  number  density  is  determined  by  assuming  a 
monodisperse  distribution  and  utilizing  Eq.  (C-2)  of  Appendix  C. 
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summarized  in  Fig.  1.  Further,  the  PDS  should  provide  real-  or  nearly  real-time  data 
acquisition  and  should  not  disturb  the  flow  field.  Hence,  the  PDS  must  provide  the  means  to 
measure  droplet  diameters  in  the  range  from  5  to  100  /tin  at  concentrations  up  to  103 
particles  per  cubic  centimeter  and  provide  continuous,  online,  in  situ  data  acquisition.  This 
set  of  criteria  was  used  to  evaluate  prospective  PDS’s  during  the  instrumentation  survey 
portion  of  the  research  program. 

Several  of  the  PDS’s  surveyed  (including  oil  slides,  rotating  multicylinders,  cascade 
impactors.  Coulter  counters,  and  commercial  optical  counters)  require  samples  to  be 
removed  from  the  flow  field.  Direct  sampling  techniques  in  which  the  particles  are  impacted 
or  otherwise  mechanically  collected  do  not,  however,  meet  the  ETF  criteria  since,  for  the 
most  part,  direct-sampling  techniques  disturb  the  flow  field  and  produce  discontinuity  in 
sampling.  Consequently,  the  majority  of  the  PDS’s  surveyed  were  electro-optical  techniques 
in  which  direct  sampling  is  not  required. 

Of  the  numerous  electro-optical  PDS’s  surveyed,  four  were  selected  to  be  evaluated 
further  in  the  icing  research  test  cell.  These  include  a  fiber-optics  sizing  system  (FOS),  two 
separate  particle-sizing  interfermometers  (PSI),  and  a  back  scattering  particle-sizing  system 
(BPS).  An  inline  holocamera  was  selected  as  the  baseline  device  for  the  experimental  part  of 
the  evaluation. 

2.2.1  Holography  System 

The  advantages  of  using  a  holocamera  as  a  particle  diagnostic  tool  are:  (1)  It  provides 
three-dimensional  viewing  capability  of  the  particle  field.  (2)  Particle  size  and  number 
density  are  obtained  directly.  (3)  When  double  pulsed,  the  holocamera  allows  obtaining 
particle  velocity  information  for  velocities  up  to  5,900  fps  (Refs.  4  through  6). 
Disadvantages  of  the  holocamera  system  are:  (1)  High  particle  number  density  can  limit 
resolution.  (2)  Data  acquisition  is  neither  continuous  nor  online.  (3)  The  data  reduction  is 
laborious,  time  consuming,  and  subjective. 

The  highest  particle  number  density  required  in  the  icing  research  cell  is  estimated  to  be 
less  than  103  cm3.  Extensive  applications  of  holography  at  AEDC  (Refs.  4  through  7)  have 
shown  that  good  quality  holograms,  as  defined  by  Royer  and  described  by  Belz  and  Menzel 
in  Ref.  7,  have  been  obtained  for  greater  particle  concentrations.  It  was  deduced,  therefore, 
that  the  particle  number  density  required  in  the  research  cell  would  not  be  large  enough  to 
severely  limit  the  resolution  of  the  holocamera.  Also,  the  time  required  to  reduce  the 
holographic  data  is  not  critical  and  does  not  severely  limit  the  holocamera’ s  application  as  a 
research  tool;  the  measurements  obtained  in  the  icing  research  test  cell  will  be  used  to 
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calibrate  the  particle  flow  field  “as  is,”  rather  than  to  set  the  flow  field  to  preselected 
conditions.  For  these  reasons  it  was  advantageous  to  use  a  holocamera  as  the  baseline 
instrument  for  AEDC  icing  research. 

Figure  4  is  a  schematic  of  the  inline  holocamera  used  in  the  icing  research  test  cell.  The 
camera  consists  of  a  double-pulsed  Apollo  ruby  laser,  a  specially  designed 
collimator/autocollimator  optical  system,  a  lens  and  filter  system,  and  a  multiple  glass  plate 
magazine  (Fig.  5).  In  addition  to  the  pulsed  laser,  a  5-mw  He-Ne  laser  was  used  to  align  the 
holography  system  before  testing.  The  holograms  were  recorded  on  Agfa-Gevaert  10E75 
high-quality  glass-backed  photographic  plates  to  minimize  image  distortion  caused  by  film 
warpage. 

The  holocamera  optical  system  is  aligned  to  image  the  centerline  of  the  icing  research  test 
cell  onto  the  film  plane  of  the  hologram.  The  sample  volume  of  the  holocamera  is  2  in.  in 
diameter  by  12  in.  in  length.  Thus,  the  holograms  contain  particle  size  information  from  a 
6-in.  range  on  either  side  of  the  film  plane.  The  holocamera  system  has  a  theoretical 
resolution  of  7.5  pm. 

Holograms  are  reconstructed  with  the  system  shown  in  Fig.  6.  The  collimated  beam  from 
a  15-mw  He-Ne  laser  provides  the  illumination  for  reconstruction.  The  hologram  is  mounted 
on  a  three-axis  traverse  system  that  allows  the  operator  to  move  the  hologram  relative  to  a 
fixed,  closed-circuit  television  vidicon.  An  auxiliary  lens  projects  an  enlargement  of  the 
reconstructed  image  in  the  vidicon  tube.  Once  an  image  is  identified  and  focused  on  the  TV 
monitor,  the  operator  indicates  the  position  of  the  image  with  a  light  pen,  and  a  Millipore® 
model  irMc  image-analysis  system  is  used  to  measure  either  the  perimeter,  diameter,  or  area 
of  the  image.  The  size  measurement  and  the  three-dimensional  position  are  recorded  on 
punched  paper  tape  for  subsequent  analysis.  Reference  7  presents  a  complete  description  of 
the  holographic  system,  including  optical  alignment  of  the  holocamera. 

2.2.2  Fiber-Optics  Particle-Sizing  System  (FOS) 

The  FOS  is  an  imaging  device  that  uses  an  expanded  laser  beam  as  a  light  source  and  an 
optical  system  to  define  a  probe  volume  in  the  particle  flow  field  (Fig.  7).  The  probe  volume 
is  focused  onto  a  linear  array  of  sensor  modules.  As  a  particle  passes  through  the  probe 
volume,  its  shadow  occludes  a  number  of  the  sensors.  The  number  of  sensor  modules 
occluded  is  proportional  to  the  particle  diameter.  Reference  8  gives  a  thorough  discussion  of 
the  theory  and  application  of  the  FOS. 

The  FOS,  as  used  in  the  icing  research  test  cell,  is  shown  schematically  in  Fig.  8.  The 
illumination  for  the  probe  volume,  located  on  the  test  cell  centerline,  is  provided  by  an  air- 
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cooled  25-mw  argon  laser.  The  laser  beam  is  spatially  filtered  and  focused  by  a  Spectra- 
Physics  collimator  and  directed  by  a  front  surface  mirror  through  a  test  cell  window. 

The  light-collecting  optics  of  the  FOS,  in  conjunction  with  a  cross-beam  coincidence 
system,  defines  the  optical  probe  volume.  Both  the  collecting  optics  and  the  coincidence 
system  are  located  outside  the  test  cell  opposite  the  argon  laser  (Fig.  8).  The  light -collecting 
optics  are  aligned  along  the  laser  axis,  and  the  cross-beam  coincidence  system  is  inclined  to 
the  laser  axis  by  12.83  deg.  The  probe  volume  is  magnified  83  times  and  projected  by  the 
collecting  optics  onto  an  array  of  32  fiber  optics,  or  “light  pipes.”  The  coincidence  system  is 
used  to  define  the  probe  volume  and  eliminate  the  spurious  effects  of  noise  coupled  with 
out-of-focus  particle  transits.  The  individual  fiber  optics,  which  are  264  jim  in  diameter,  are 
arrayed  along  a  straight  line  for  a  total  length  of  0.3325  in.;  each  fiber  terminates  at  a 
photomultiplier  (PM)  tube.  Each  of  the  32  photodetector  circuits  provides  information  to 
the  size-determination  logic  of  the  data  handling  system,  as  described  in  Ref.  9.  The  number 
of  events  of  specified  size  is  stored  and  displayed  periodically  on  both  a  cathode  ray  tube 
(CRT)  and  a  digital  display  on  the  front  panel  of  the  particle  data  system  (Fig.  9).  A 
Berkely®  EPUT  meter  is  slaved  to  the  memory-enable-switch  of  the  data  acquisition  system 
so  that  the  time  over  which  the  measurements  are  made  can  be  accurately  recorded. 

Thirty  channels  of  particle  size  data  are  stored  by  the  data  handling  system.  Data  output 
is  controlled  manually.  The  channel  number  and  the  particle  counts  corresponding  to  that 
channel  comprise  the  output  data;  each  is  primed  out  sequentially  by  a  Hewlett  Packard 
5055A  digital  recorder. 

A  complete  description  of  the  FOS  as  used  in  the  icing  research  test  cell  is  given  in  Ref.  9. 

2,2.3  Particle-Sizing  Interferometer  (PSI) 

The  principle  of  operation  of  the  PSI  is  illustrated  in  Fig.  10.  Two  laser  beams  of  equal 
intensity  and  linearly  polarized  in  the  same  direction  are  focused  and  crossed  to  form  a  small 
ellipsoidal  probe  volume  in  space.  Because  of  the  interaction  of  the  two  beams,  the  probe 
volume  is  filled  with  parallel,  equally  spaced  fringes.  As  a  particle  passes  through  the  probe 
volume,  light  is  scattered  which  produces,  by  means  of  a  photomultiplier  tube,  a  current 
signature  characteristic.  As  shown  in  Fig.  10,  the  current  waveform  is  composed  of  a  high- 
frequency  A-C  component  and  a  Gaussian-shaped  “pedestal”  component.  The  ratio  of  A-C 
amplitude  divided  by  the  pedestal  amplitude  is  defined  as  the  visibility  of  the  waveform  and, 
as  shown  in  Ref.  10,  the  particle  size  is  a  unique  factor  (within  limits)  of  the  visibility. 

Hence,  by  utilizing  the  signature  characteristics,  the  size  of  each  particle  that  passes 
through  the  probe  volume  can  be  determined.  More  detailed  descriptions  of  interferometric 
particle  sizing  are  given  in  Refs.  10  and  11. 
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The  major  subsystems  of  a  PSI  may  be  categorized  as  (1)  a  laser  (which  provides  a 
coherent,  highly  collimated  light  source),  (2)  transmitting  optics  (which  precondition  and 
focus  the  laser  light  to  form  the  probe  volume),  (3)  receiving  optics  (which  collect  and  focus 
the  light  scattered  by  particles  moving  through  the  probe  volume),  and  (4)  signal 
conditioning  electronics  (which  amplify  the  scattered  radiation  and  process  the  signal  to 
obtain  visibility).  The  specification  and  arrangement  of  these  subsystems’  components  are 
dictated  by  the  intended  application  of  the  PSI. 

Two  PSI  systems  were  selected  for  evaluation  in  the  icing  research  test  cell.  One  system  is 
a  PSI  developed  at  AEDC  and  the  other  is  a  commercially  available  system.  Although  the 
operating  principle  is  the  same  for  both  systems,  the  differences  existing  between  them  are 
great  enough  that  experimental  evaluation  of  each  was  required.  A  description  of  each  PSI 
as  tested  in  the  icing  research  test  cell  is  given  below. 

2.2.3. 1  AEDC-Developed  PSI 

The  optical  system  is  shown  schematically  in  Fig.  1 1  (Fig.  1 1  and  much  of  the  following 
descriptions  were  taken  from  Ref.  12).  The  light  source  is  a  15-Mw  helium-neon  laser 
operating  in  the  TEMoo  mode  at  a  wavelength  of  0.6328  ftm.  The  transmitting  optics 
consist  of  a  Bragg  cell  beamsplitter  driven  at  15  MHz  and  a  system  of  lenses  which  brings 
the  two  beams  to  simultaneous  cross  focus;  the  crosspoint  is  found  39.4  in.  in  front  of  the 
final  transmitting  lens.  Use  of  a  Bragg  cell  as  a  beamsplitter  produces  a  set  of  moving  fringes 
in  the  probe  volume,  and  the  Bragg  cell  is  driven  in  a  manner  to  move  the  fringes  in  the 
direction  opposite  the  direction  of  the  droplet  flow.  A  fringe  spacing  of  19  gm  was  used; 
however,  other  fringe  spaces  can  be  selected  by  changing  the  transmitting  lenses.  The 
beamwaist  diameter  is  150 /*m.  A  set  of  5.9-in.-diam,  f/5  lenses  is  used  to  collect  and  focus 
the  scattered  light.  A  special  mask  that  admits  light  through  two  60-deg  wedge-shaped 
openings  is  used  on  the  collection  lens  (see  Fig.  1 1).  A  200-jtm  pinhole  in  front  of  the  PM 
tube  blocks  light  originating  at  points  other  than  the  probe  volume.  The  current  waveforms 
are  produced  with  an  RCA®  8644  PM  tube. 

The  electronics  used  to  measure  the  visibility  are  shown  as  a  block  diagram  in  Fig.  12. 
The  photodetector  signal  is  amplified  and  sent  simultaneously  to  a  bandpass  filter  and  to  a 
low-pass  filter  for  separation  of  the  A-C  component  from  the  pedestal.  Then,  the  pedestal 
and  A-C  magnitude  are  peak  detected.  The  peak  pedestal  signal  is  fed  to  the  denominator 
input  of  analog  divider  while  the  peak  of  the  A-C  signal  fed  to  the  divider  is  numerator  input 
of  the  divider.  (The  output  signal  from  the  divider  is  ratio  of  peak  A-C  magnitude  to  peak 
pedestal.) 

The  data  handling  system  consists  of  a  microcomputer  that  stores  the  data  in  histogram 
form  and  a  Tektronix  Model  4051  table-top  computer  that  controls  data  acquisition.  It  has 
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been  noted  (Ref.  12)  that  the  particle  size  data  from  the  PS1  should  be  corrected  to 
compensate  for  the  size  of  the  sample  volume,  which  varies  with  particle  size,  and  also  to 
account  for  variations  of  the  bin  width  of  the  histograms,  which  vary  with  visibility. 
Consequently,  the  40S1  is  also  used  to  correct  the  raw  data  from  these  variations  so  that  the 
bar  heights  of  the  histograms  are  output  as  number  of  particles  per  cubic  centimeter  per 
micrometer.  The  output  data  are  comprised  of  histogram  bar  heights  as  a  function  of 
particle  diameter,  with  a  maximum  resolution  of  32  bins  per  histogram. 

2.2.3.2  Commercial-PSI 

The  commercially  available  PSI  was  built  to  AEDC  specifications  by  Spectron 
Development  Laboratories,  Inc.  The  optical  system  is  shown  in  Fig.  13.  The  light  source  is  a 
15-mw  He-Ne  laser  operating  in  the  TEMoo  mode  at  a  wavelength  of  0.6328  /im.  The 
transmitting  optics  consist  of  a  glass  block  beamsplitter  and  a  lens  system  that  cross  and 
focus  the  beams  at  the  same  location.  The  crosspoint  of  the  beams  is  41  in.  in  front  of  the 
lens.  A  fringe  spacing  of  50  **m  was  used  during  the  experimental  investigation.  The  beam- 
waist  diameter  is  500  fii n.  The  collection  optics  consist  of  a  set  of  5-in.-diam,  f/5  lenses;  a 
90-deg  mask  (Fig.  13)  is  used  on  the  collection  lens.  An  EMI9781R  side  window  PM  tube  is 
used  as  the  light  detector. 

The  PSI  electronics  consist  of  a  preamplifier  contained  with  the  photomultiplier,  a 
Krohn-Hite  model  3202  variable  filter,  and  a  PMAP  550  visibility/velocity  processor  (VVP) 
(Ref.  13),  The  electronics  are  shown  in  block  diagram  form  in  Fig.  14.  Signals  from  the  PMT 
are  divided  into  two  portions,  the  first  of  which  is  sent  to  an  integrator  and  integrated 
directly.  The  integral  -  called  the  pedestal  integral  -  yields  a  number  which  is  directly 
proportional  to  the  mean  value  of  the  amplitude  since  the  sinusoidal  component  of  the 
signal  will  average  to  zero.  The  value  of  the  integral  is  then  fed  to  the  denominator  input  of 
an  analog  divider.  The  second  portion  of  the  signal  is  bandpass  filtered  to  remove  the 
pedestal  and  the  filtered  output  is  full  wave  rectified  then  integrated.  The  integrated  value  of 
the  A-C  signal  is  fed  into  the  numerator  input  of  the  divider.  The  output  from  the  divider  is 
sent  to  an  analog-to-digital  converter. 

The  data  handling  system  is  comprised  of  a  microprocessor,  high-resolution  CRT 
monitor,  and  a  printer  with  both  alphanumerics  and  graphics  capabilities.  The  system 
displays  the  particle  size  data  in  histogram  form  at  a  maximum  resolution  of  64  bins  per 
histogram.  The  particle  size  data  are  reduced  to  prescribed  units  before  printout;  for  the 
present  investigation,  the  output  was  presented  in  number  of  particles  as  a  function  of 
particle  diameter. 
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2.2.4  Backscattering  Particle-Sizing  System  (BPS) 

When  a  particle  is  illuminated,  it  interacts  with  the  electromagnetic  radiation  and  scatters 
light  in  all  directions.  The  intensity  of  the  scattered  light  is  directly  related  to  the  intensity  of 
the  original  illumination,  the  effective  size  of  the  particle,  and  the  surface  characteristics  of 
the  particle.  Several  instruments  have  been  built  which  rely  on  the  principle  of  light 
scattering  to  size  particles  one  at  a  time.  Such  a  device  is  shown  in  Fig.  15.  A  particle  in  the 
sampling  volume  scatters  light  which  is  collected,  detected,  and  converted  into  a  voltage. 
The  voltage,  above  a  residual  background  value,  is  related  through  a  calibration  to  the  size 
of  the  particle.  Size  distribution  is  determined  by  counting  the  number  of  pulses  of  given 
magnitudes.  If  the  pulse  rate  and  the  particle  field  velocity  are  known,  then  the  particle 
concentration  can  be  determined. 

The  light  scattering  device  selected  for  evaluation  in  the  icing  research  test  cell  was  built 
by  Environmental  Systems  Corporation.  The  device  uses  backscattered  light  to  size  particles; 
that  is,  in  Fig.  1 5  the  angle  6  between  the  laser  beam  and  the  light  collection  optics  is  1 80  deg. 
The  beam  forming  and  light  collecting  optics  are  arranged  co-axially,  and  are  located  in  a 
3- 1/2-in. -diam  tube.  The  tube  protrudes  into  and  disturbs  the  flow  Field  (Fig.  16).  Since  it 
disturbs  the  flow  field,  the  BPS  does  not  meet  all  of  the  requirements  of  the  ETF  (see  Section 
2.2.0).  However,  it  was  felt  that  experimental  evaluation  would  provide  valuable  insight  into 
the  particle-sizing  capabilities  of  such  a  device. 

The  optical  system  is  shown  in  Fig.  17.  The  light  source  is  a  Ga-As  PN-junction  laser  that 
emits  a  200-nsec-wide  light  pulse  at  103  pulses/sec.  The  wavelength  of  the  laser  light  is  0.904 
pm.  The  emitted  light  is  focused  by  a  three-lens  system  to  form  the  sampling  volume.  The 
sample  volume  is  0.5  mm3  and  is  located  6  in.  in  front  of  the  lens  system.  The  light  scattered 
by  particles  in  the  sample  volume  is  collected  by  a  pair  of  1-1/2-in. -diam  f/6  lenses.  The 
scattered  light  is  divided  into  equal  parts  by  a  mirror-type  beamsplitter  and  then  focused 
onto  two  separate  detector-preamplifier  assemblies.  The  two  detectors  along  with  proper 
apertures  provide  a  particle  rejection  system  that  ensures  that  a  particle  passes  through  the 
calibration  portion  of  the  sample  volume. 

The  signals  from  the  detectors  are  amplified,  digitized,  and  processed  in  a 
microprocessor.  The  signals  from  the  two  detectors  are  first  compared  to  ensure  that  they 
originate  from  a  particle  inside  the  specified  portion  of  the  sample  volume.  Qualified  signals 
are  then  classified  according  to  magnitude  and  counted  in  eight  preset  particle  size  bins. 
Data  output  is  manually  controlled.  The  channel  number  and  the  particle  counts 
corresponding  to  that  channel  comprise  the  output  data.  The  particle  count  is  manually 
recorded  from  the  LED  display  located  on  the  control  panel  (Fig.  16). 


13 


AEDC-TR  -80-33 


2.3  SUPPORTING  INSTRUMENTATION 

Conventional  instrumentation  was  required  to  determine  test  conditions  during  the 
experimental  program.  The  instrumentation  stations  are  indicated  in  Fig.  3.  Total 
temperatures  were  measured  with  single-shielded,  self-aspirating,  copper-constantan 
thermocouple  probes  with  an  ambient  reference  junction.  Pressures  were  measured  with 
strain-gage-type  absolute  pressure  transducers  referenced  to  atmospheric  pressure. 

» 

A  combination  total/static  pressure  probe  (Fig.  3a,  detail  A)  was  installed  in  the  test 
section.  Radial  surveys  of  total/static  pressure  were  made  across  the  test  section  with  the 
probe.  The  probe  positions  were  determined  with  an  absolute  BCD  encoder  connected  to 
the  probe. 

A  turbine-type  flowmeter,  installed  in  the  water  line  to  the  spray  nozzle,  was  used  to 
determine  the  water  content  of  the  icing  cloud.  This  flowmeter  measured  all  water  flowing 
into  the  test  cell  through  the  spray  nozzle  and  provided  primary  water  content  data.  The 
turbine  meter  has  a  dynamic  range  of  from  2,5  to  21  gallons  per  hour. 

A  rotameter,  used  to  determine  the  airflow  rate  to  the  spray  nozzle,  was  installed  in  the 
atomizing  air  line  to  the  spray  nozzle.  The  rotameter  is  rated  at  300  psi  and  has  a  dynamic 
range  for  airflows  of  from  0.5  to  10  standard  cubic  feet  per  min. 

2.4  TEST  PROCEDURES 

All  instrumentation  parameters  of  the  support  instrument  were  calibrated  before  each 
test.  The  calibrations  were  conducted  in  place  with  ambient  pressure  in  the  test  cell.  The 
calibrations  included  millivolt  calibration  of  all  temperature  parameters  and  resistance 
calibration  of  all  pressure  parameters.  All  calibrations  were  reviewed  and  any  discrepancies 
were  corrected  before  the  test. 

After  calibrations  were  performed,  the  test  cell  pressure  was  reduced  to  the  desired 
pressure  altitude.  Cell  conditions  were  allowed  to  stabilize  for  approximately  five  minutes, 
and  a  data  point  was  obtained  for  final  instrumentation  check.  During  cell  stabilization,  a 
final  alignment  check  was  made  on  the  particle  diagnostic  system  to  be  tested. 

After  the  final  instrumentation  check,  both  primary  and  secondary  inlet  air  was  admitted 
to  the  test  cell  at  the  required  pressure  and  temperature.  Airflow  through  the  spray  nozzle 
was  also  started  and  was  maintained  throughout  the  test  to  prevent  freezing  of  the  spray 
nozzles.  Once  the  test  cell  flow  conditions  had  stabilized  to  steady  state,  the  air  purge 
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through  the  water  outlet  of  the  spray  nozzle  was  stopped,  and  the  desired  water  and  air  flows 
through  the  spray  nozzle  were  set.  The  spray  water  flow  rates  were  selected  so  that  for  a 
given  bulk  velocity,  the  resulting  liquid  water  content  (LWQ  varied  with  the  meterorological 
conditions  of  Fig.  1  (i.e.,  LWC  varied  from  0.2  to  3  gm/m3).  Droplet  size  information  was 
then  obtained  with  use  of  the  PDS,  and  at  least  two  holograms  or  histograms  were  taken  for 
each  spray  nozzle  test  condition.  A  high-speed  minicomputer  was  used  to  record  the  test  cell 
data  reported  concurrently  with  the  hologram/histogram  data.  When  the  test  sequence 
ended,  the  PDS  was  reset,  the  next  spray  nozzle  condition  established,  and  the  data 
acquisition  process  repeated. 

The  matrix  of  test  conditions  utilized  in  the  experimental  program  is  presented  in  Table 
2.  The  matrix  was  chosen  to  provide  the  largest  range  of  particle  sizes  over  which  the  PDS 
must  operate  during  icing  simulation. 

The  test  procedure  was  utilized  throughout  the  experimental  program  to  evaluate  each 
PDS.  The  holocamera  was  first  used  to  determine  the  water  droplet  size  distribution  for  each 
test  condition  summarized  in  Table  2.  Each  of  the  remaining  PDS’s  was  then  tested  to 
obtain  water  droplet  size  distributions  for  each  test  condition  of  the  experimental  program. 

2.5  PRECISION  OF  MEASUREMENTS 

Uncertainties  (bands  which  include  95  percent  of  the  calibration  data)  of  the  basic 
experimental  parameters  were  estimated  from  repeat  calibrations  of  the  instrumentation 
systems.  These  uncertainties  were  estimated  against  secondary  standards  whose  precisions 
were  traceable  to  the  NBS.  The  uncertainties  were  combined  by  means  of  the  Taylor  series 
method  of  error  propagation  (Ref.  14)  to  determine  the  precision  of  the  experimental 
parameters  presented  in  Table  3. 

3.0  RESULTS  AND  DISCUSSION 

3.1  HOLOGRAPHIC  DATA 

3.1.1  Holographic  System  Calibration 

Calibrations  of  the  holographic  system  (i.e.,  the  holocamera  and  reconstruction  system) 
were  conducted  periodically  throughout  the  test  program.  The  purposes  of  the  calibrations 
were  to  determine  the  diameter  of  the  smallest  droplet  resolvable  with  the  holographic 
system  and  to  establish  a  set  of  calibration  factors  for  the  system.  Repeat  calibrations  of  the 
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system  were  necessary  throughout  the  program  to  ensure  that  the  correct  calibration  factors 
were  used  to  reduce  the  holographic  data. 

An  Ealing  resolution  chart  was  used  to  accomplish  the  calibration.  The  resolution  chart 
is  a  dot  and  line  comparison  reticle  that  is  made  up  of  two  rows  of  dots  and  lines  of  equal 
width;  the  sizes  vary  from  3  pm  to  100  pm.  The  resolution  chart  was  placed  at  various  radial 
locations  in  the  test  plane  of  the  research  test  cell.  Radial  locations  of  1,6,  7.5,  14.5,  and 
17.5  cm  on  either  side  of  the  test  cell  centerline  were  chosen  to  ensure  that  calibrations  were 
obtained  for  the  entire  droplet  flow  field.  Two  holograms  were  taken  of  the  resolution  chart 
at  each  radial  location  for  a  total  of  10  holograms  per  calibration.  Each  hologram  was 
reconstructed  to  obtain  the  diameter  size  of  each  of  the  removable  dots. 

By  plotting  the  actual  dot  diameters  as  functions  of  the  measured  dot  diameters,  the 
systems  calibration  function  was  obtained;  this  is  illustrated  in  Fig.  18.  Each  value  of 
measured  diameter  shown  in  the  figure  represents  an  average  of  10  diameters.  The  error  bars 
shown  in  Fig.  18  represent  the  highest  and  lowest  values  of  measured  diameter  taken  from 
the  10  sets  of  data. 

Figure  18  shows  that  the  smallest  droplet  resolvable  with  the  holographic  system  is  10 
/an.  Also,  the  calibration  function,  represented  by  the  straight  line  fit  to  the  calibration  data, 
apparently  changed  during  the  experimental  program. 

The  calibration  functions  were  determined  from  a  linear  regression  curve  Fit  to  the 
calibration  data.  In  all  three  cases,  the  coefficient  of  determination,  r2,  indicates  that  the 
quality  of  Fit  achieved  by  the  regression  was  0.997  or  better.  (Note  that  a  value  close  to  1.00 
indicates  a  better  fit  than  values  close  to  zero.)  The  equations  used  to  calculate  the  regression 
coefficients  and  the  coefficient  of  determination  are  given  in  Appendix  B. 

The  largest  shift  in  calibration  occurred  during  the  First  portion  of  the  experimental 
program,  between  the  pretest  calibration  and  the  calibration  during  testing.  A  comparison 
of  these  calibrations  revealed  that  the  absolute  error  in  droplet  diameter  as  a  result  of  change 
in  calibrations  varied  from  1 .83  /mi  at  the  smallest  (10  pm)  droplet  diameter  to  5.26 pm  at  the 
largest  (100  pm)  diameter.  Hence,  use  of  the  incorrect  calibration  function  to  reduce  the 
holographic  data  during  this  portion  of  the  experimental  program  would  increase  the  value 
of  the  smallest  resolvable  droplet  diameter  from  10  pm  to  1 1 .83  pm,  and  in  a  like  manner  the 
value  of  the  largest  droplet  diameter  would  be  increased  from  100  pm  to  105.26  pm.  Care 
was  therefore  taken  throughout  the  experimental  program  to  ensure  that  the  correct 
calibration  was  used  to  reduce  the  holographic  data. 


16 


AEDC-TR-80-33 


3.1.2  Research  Cell  Calibration 

The  icing  research  test  cell  was  calibrated  with  the  holographic  system.  With  the  research 
test  cell  in  the  basic  configuration  (16  ft  axial  distance  between  the  bellmouth  entrance  and 
the  test  plane)  (Fig.  3A),  holograms  were  taken  at  each  of  the  spray  nozzle  conditions  listed 
in  Table  2.  From  the  28  holograms  taken  of  the  14  spray  nozzle  conditions,  7,250  droplet 
diameters  were  measured.  From  200  to  300  droplets  were  measured  from  each  hologram.  It 
was  assumed  that  400  droplet  counts  were  sufficient  for  statistical  analysis  of  each  spray 
nozzle  condition. 

The  droplet  size  information  obtained  from  the  hologram  was  used  to°  construct 
histograms  of  droplet  size  distribution.  The  histograms  are  plots  of  the  frequency  of 
occurrence  as  a  function  of  the  droplet  size  range.  The  histograms  shown  in  Figs.  19, 20,  and 
21  typify  the  calibration  data. 

The  calibration  data  were  also  used  to  construct  accumulative  distributions  of  the 
percent  of  droplets  greater  than  a  given  droplet  size;  the  data  were  plotted  against  droplet 
size  on  log-probability  paper.  The  log-probability  plots  obtained  for  the  histograms  of  Figs. 
19  through  21  are  shown  in  Figs.  22  through  24,  respectively.  Examination  of  Figs.  22 
through  24  reveals  that  the  plotted  data  lie  on  an  essentially  straight  line.  This  indicates  that 
the  calibration  data  are  logarithmically  distributed  (Ref.  14).  Therefore,  a  normal 
probability  equation  with  a  logarithmic  variate  can  be  applied  to  the  droplet  size 
distributions  of  the  calibration  data.  The  frequency  with  which  a  droplet  of  a  given  diameter 
occurs  is 
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The  frequency  distribution  is  a  unique  function  of  the  geometric  mean  size,  M,  and  the 
geometric  standard  deviation,  Sg;  once  they  have  been  specified,  only  one  distribution  is 
possible.  To  compare  the  droplet  distribution  obtained  from  calibration  with  those  droplet 
distributions  obtained  from  calibration  with  those  droplet  distributions  obtained  by  other 
particle-diagnostic  systems,  only  the  values  of  M  and  Se  for  both  distributions  need  be 
compared.  The  geometric  mean  size  is  given  by 


fog  M 
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and  the  geometric  standard  deviation  by 
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Although  the  values  of  M  and  Sg  can  be  calculated  directly  from  the  histogram  data  with 
Eqs.  (2)  and  (3),  these  measures  can  easily  be  obtained  graphically  from  the  cumulative 
probability  plots.  When  droplet  diameters  are  log-normally  distributed,  the  median  droplet 
size  is  equal  to  the  geometric  mean  droplet  size  (Ref.  14).  The  value  of  M  is,  for  a  given 
distribution,  the  droplet  diameter  corresponding  to  the  median  (50-percent  value)  of  the 
probability  coordinate  of  the  cumulative  probability  plots.  Further,  it  can  be  shown  (Ref. 
14),  that  the  geometric  standard  deviation  can  be  obtained  from  the  plots  as  the  ratio  of  the 
size  at  15.87-percent  probability  to  that  at  50  percent,  or  the  ratio  at  50-percent  probability 
to  that  at  84.13  percent.  Both  the  analytical  method  utilizing  Eqs.  (2)  and  (3)  and  the 
graphical  method  just  described  were  used  to  obtain  values  of  the  geometric  mean  diameter 
and  geometric  standard  deviation.  Table  4  presents,  for  comparison,  values  of  M  and  Sg 
determined  by  both  methods.  The  geometric  mean  size  defined  by  Eq.  (2)  or  as  determined 
with  the  graphical  method  is  based  on  the  count  frequency  of  the  given  diameter. 
Historically  the  mass  median  diameter  (dgm)  is  the  central  tendency  specified  as  the  “mean 
effective”  droplet  diameter  of  an  icing  cloud  (Refs.  15  and  16).  This  diameter  is  a  mass- 
weighted  average  and  is  defined  as  the  diameter  of  those  droplets  for  which  one-half  of  the 
entire  liquid  water  content  of  the  icing  cloud  is  contained  in  droplets  whose  diameters  are 
less  than  dBm.  Conversion  of  the  geometric  mean  diameter,  M,  into  the  mass  median 
diameter,  dgm.  is  accomplished  by  the  transformation 

f!ndgm  =  ?n  iVl  +  3<('n  S^2  W 


Use  of  Eq.  (4)  is  appropriate  only  if  a  log-normal  distribution  exists  (Ref.  14).  Since  it  had 
been,  shown  that  all  of  the  calibration  data  are  indeed  log-normally  distributed,  then  Eq.  (4) 
was  used  to  determine  the  mass  median  diameters  reported  in  Table  4.  The  standard 
deviations  in  Table  4  are  the  geometric  standard  deviations.  As  pointed  out  in  Ref.  14,  one 
favorable  aspect  of  the  log-normal  distribution  is  that  regardless  of  the  basis  used  for  its 
determination,  the  geometric  standard  deviation  is  the  same.  Consequently,  the  results  listed 
in  Table  4  completely  define  the  holographic  calibration  of  the  spray  nozzle  conditions  of 
Table  2,  for  comparison  either  on  a  number  basis  or  on  a  mass  basis. 
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The  mass  median  diameters  listed  in  Table  4  were  used  in  constructing  a  calibration  curve 
for  the  spray  nozzle.  The  calibration  curve  (Fig.  25)  is  a  plot  of  mass  median  diameter  as  a 
function  of  the  spray  nozzle  mass  flow  ratio.  (The  mass  flow  ratio  is  defined  as  the  mass 
flow  rate  of  spray  nozzle  atomizing  air  divided  by  the  mass  flow  rate  of  spray  nozzle  water.) 
A  power  curve  of  the  form  y  =  axk  was  fitted  to  the  holographic  data  and  is  shown  as  the 
solid  line  in  Fig.  25.  The  regression  coefficients  for  the  curve  fit  were  determined  by  the 
least-squares  method.  (The  equations  used  are  given  in  Appendix  B.)  The  curve-fit  equation 
is 


d 


gm 


17.87  (*./*, 


-0.3  IB 


(5) 


The  coefficient  of  determination  for  the  fit  is  0.983,  which  indicates  an  excellent  fit  to  the 
data.  The  dashed  lines  of  Fig.  25  serve  to  illustrate  that  all  of  the  calibration  data  fall  within 
±■2  pm  of  the  curve. 

The  spray  nozzle  operating  characteristics  (Fig.  25)  are  fairly  “flat”  in  the  15-  to  25-pm 
range;  large  changes  in  the  mass  flow  ratio  are  required  to  produce  small  changes  in  the  mass 
median  droplet  diameters.  Hence,  in  this  range  the  measured  flow  rates  of  spray  nozzle 
atomizing  air  and  water  can  be  used  with  a  high  level  of  confidence  to  determine  the  mass 
median  diameter  of  the  icing  clouds. 

This  is  not  the  case,  however,  for  droplet  diameters  larger  than  approximately  30  pm,  for 
in  this  region  of  the  curve,  small  changes  in  the  mass  flow  ratio  produce  large  changes  in  the 
mass  median  diameter.  Hence,  when  the  spray  nozzle  is  used  to  produce  a  mass  median 
droplet  diameter  larger  than  30  pm,  the  calibration  curve  shown  in  Fig.  25  cannot  be  used; 
rather  an  independent  droplet  measurement  system,  such  as  a  holocamera  or  one  of  the 
other  PDS’s  described  in  Section  2.0,  is  required  to  ensure  accurate  control  of  dgn,t 

The  holographic  calibration  data  were  also  used  to  determine  the  LWC  at  the  test  section 
of  the  research  test  cell.  The  LWC  was  calculated  by  dividing  the  weight  of  all  droplets 
counted  by  the  volume  of  the  test  section  scanned  during  data  reduction;  the  equation  is 
given  in  Appendix  A.  In  Fig.  26  the  LWC  values  measured  at  the  spray  nozzle  are  compared 
with  the  values  calculated  from  the  holographic  data.  The  data  clearly  demonstrate  that  the 
calculated  LWC  is  biased  toward  higher  values.  Two  conditions  that  contribute  to  the 
biasing  are  (1)  improper  droplet  focus  recognition  during  hologram  reconstruction  and  (2) 
large  snow  or  ice  crystals  in  the  flow  field.  Since  LWC  depends  on  d?,  the  error  in  the 
calculated  LWC  will  be  about  three  times  larger  than  the  corresponding  error  made  during 
measurement  of  the  droplet  diameter;  hence,  failure  to  correctly  size  droplets  causes  serious 
errors  in  the  subsequent  LWC  calculations. 
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Droplet  focus  recognition  is  based  upon  the  rate  of  change  of  light  intensity  near  the 
droplet  edges  (Fig.  27).  If  the  edges  are  not  in  focus,  the  droplet  will  appear  larger  than  it 
should,  biasing  the  droplet  sizes  toward  higher  values.  When  snow  and/or  ice  crystals  are 
present  in  the  flow  field,  the  density  of  water  is  used  to  calculate  the  droplet  weight  rather 
than  the  snow/ice  density,  and  the  assumption  of  spherical  droplets  inherent  to  the  LWC 
calculations  may  be  violated,  similarly  biasing  the  LWC  toward  higher  values.  Determining 
LWC  from  holographic  data  is  not,  at  this  time,  a  viable  means  of  obtaining  an  accurate 
measure. 

3,2  PARTICLE  DIAGNOSTIC  SYSTEM  DATA 

Droplet  size  measurements  are  presented  for  each  PDS  tested.  Data  were  obtained  for 
each  PDS  at  the  test  conditions  listed  in  Table  2.  The  research  test  cell  was  used  in  the  basic 
configuration,  that  is,  with  a  16-ft  axial  distance  from  the  bellmouth  entrance  to  the 
measurement  plane  of  the  PDS.  At  least  two  droplet  size  data  sets  were  obtained  with  each 
PDS  for  each  test  condition  (Table  2),  and  at  least  1,000  droplets  were  sized  per  data  set.  In 
several  cases,  multiple  data  sets  were  made  at  given  spray  nozzle  test  conditions  to  verify 
repeatability  of  the  results  or  to  correct  experimental  deficiencies. 

3.2.1  FOS  Results 

Figure  28  is  a  typical  histogram  constructed  from  the  droplet  size  data  obtained  with  the 
FOS.  The  first  two  bins,  centered  at  droplet  diameter  3.065  /an  and  6.13  #tm,  respectively, 
are  empty.  The  data  channels  corresponding  to  these  two  bins  were  extremely  noisy,  and 
spurious  counts  were  recorded  in  these  bins  during  data  acquisition.  The  noise  was 
determined  to  be  electro-optically  induced  (i.e.,  caused  by  weak  or  out-of-focus  droplet 
images).  Since  attempts  to  eliminate  the  noise  in  channels  1  and  2  were  unsuccessful,  these 
two  channels  were  “amputated”  from  the  histograms.  Hence,  the  smallest  droplet 
resolvable  with  the  FOS,  as  used  in  this  research,  is  approximately  7.5  /tm — the  smallest 
droplet  diameter  corresponding  to  the  third  bin. 

Figure  29  is  a  log-probability  plot  of  the  histogram  data  shown  in  Fig.  28.  The  data  lie 
along  a  straight  line,  indicating  that  the  droplet  diameters  are  log-normally  distributed.  This 
was  the  case  for  all  droplet  size  data  taken  with  the  FOS. 

The  distribution  of  the  mass  median  droplet  diameters  obtained  from  the  FOS  data  is 
shown  in  Fig.  30  as  a  function  of  the  spray  nozzle  mass  ratio.  The  results  are  compared  to 
the  empirical  distribution  for  the  mass  median  diameters  obtained  from  the  holographic 
data  [Eq.  (5),  Section  3.12].  This  expression  provides  a  reasonable  fit  to  the  data;  most  of 
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the  FOS  results  lie  within  the  error  band  of  the  holographic  data,  the  dashed  lines  of  Fig.  30. 
This  error  band  is  plus  and  minus  two  standard  deviations  { ±  2  Sg)  of  the  variation  between 
the  holocamera  data  and  the  empirical  distribution  of  Eq.  (5).  (See  Appendix  B  for  the 
determination  of  Sg.) 

To  determine  the  accuracy  with  which  Eq.  (5)  predicts  the  FOS  data,  the  deviations 
between  the  measured  droplet  diameters  and  the  predicted  diameters  from  Eq.  (5)  are 
plotted  as  a  function  of  droplet  diameter  (Fig.  31).  In  Fig.  31  the  deviations  have  been  fitted 
with  a  linear  regression  least-squares  fit.  The  least-squares  line  deviates  slightly  from  the  line 
of  zero  error.  However,  over  the  region  of  interest  (dgm  =  15  pm  to  dgn)  =  25  nm)  the 
deviation  is  less  than  one  micrometer,  which  is  well  within  the  data  scatter  of  the 
holocamera.  Therefore,  the  FOS  can  size  icing  cloud  droplets  as  accurately  as  the 
holocamera. 

3.2.2  AEDC-PS1  Results 

A  typical  histogram  of  the  droplet  size  information  obtained  with  the  AEDC-PSI  is 
shown  in  Fig.  32.  This  histogram  differs  in  shape  from  the  histograms  obtained  with  the 
holocamers  data  (Figs.  19  through  21)  in  that  at  least  two  maxima  or  modes  of  the 
frequency,  ft,  are  exhibited.  Two  modes  correspond  to  droplet  diameters  of  12.5  gm  and  22 
pm.  Figure  33  is  a  log  probability  plot  of  this  typical  histogram.  The  data  do  not  plot  as  a 
straight  line;  hence,  unlike  the  holocamera  data,  the  droplet  size  distributions  obtained  with 
the  AEDC-PSI  are  not  log-normal. 

Irani  and  Callis  (Ref.  17)  note  that  the  distribution  shown  in  Fig.  33  indicates  a  log¬ 
normal  distribution  from  which  particles  have  been  removed.  To  determine  whether  the 
apparent  loss  in  droplets  was  caused  by  changes  in  the  research  test  cell  operating 
characteristic  or  whether  the  AEDC-PSI  incorrectly  sized  the  droplets,  PSI  data  taken  on 
different  days  for  the  same  spray  nozzle  test  conditions  were  compared.  Figure  34  is  a  typical 
log  probability  plot  showing  the  comparison.  As  seen  in  Fig.  34,  the  results  obtained  with  the 
PSI  correlate  well  for  the  two  test  dates  and  are  self  consistent. 

Nine  of  the  data  points  reported  for  the  holocamera  in  Sec.  3.1  and  all  of  the  FOS  data 
reported  in  Sec.  3.2.1  were  taken  in  the  interval  between  the  PSI  test  dates,  and  all  of  the 
distributions  obtained  with  both  devices  were  found  to  be  log-normal.  Because  of  the 
excellent  correlation  between  the  holocamera  data  and  FOS  data,  it  was  concluded  that  the 
droplet  clouds  produced  in  the  research  test  cell  are  log-normally  distributed.  Further, 
because  the  data  obtained  with  the  AEDC-PSI  were  found  to  exhibit  abnormal  log-normal 
yet  self-consistent  distributions,  it  was  concluded  that  the  apparent  loss  of  droplets  indicated 
by  the  log-probability  plots  was  attributable  to  measurement  inaccuracies  of  the  AEDC-PSI. 
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The  inaccuracies  in  the  droplet  size  data  produced  by  the  AEDC-PSI  apparently  result 
from  the  relatively  high  number  density  that  exists  in  a  simulated  icing  cloud.  The  existence 
of  many  droplets  in  the  measurement  plane  increases  the  probability  that  multi-droplets  are 
present  in  the  sample  volume  during  data  acquisition.  Also,  the  background  noise  which  is 
produced  by  droplets  that  cross  the  laser  beams  outside  of  the  crossover  region  is  increased. 

Either  condition  or  both  would  result  in  biased  droplet  size  data  (Ref.  18). 

* 

In  the  case  of  the  AEDC-PSI,  the  data  bias  produced  by  large  number  densities  is  toward 
large  droplet  diameters.  This  is  illustrated  in  Fig.  35,  which  is  a  plot  of  the  geometric  mass 
mean  droplet  diameters  as  a  function  of  the  spray  nozzle  mass  ratio.  Equation  (5)  is  also 
plotted  in  Fig.  35  for  comparison.  As  seen  in  Fig.  35  the  dgir,  values  for  the  AEDC-PSI  data 
are  insensitive  to  variation  in  spray  nozzle  mass  ratio;  also,  in  the  range  of  interest  (15  to 
25  /tm),  they  exceed  the  dgm  predicted  by  the  holocamera  data.  It  is  concluded,  therefore, 
that  the  AEDC-PSI  cannot  be  used,  as  presently  designed,  to  size  droplets  in  a  simulated 
icing  cloud. 

3.2.3  Commercial-PSl  Results 

Figure  36  is  a  typical  histogram  of  the  droplet  size  data  obtained  with  the  commercial- 
PSI.  Compared  to  the  holocamera  results  shown  in  Fig.  19,  the  histogram  of  Fig.  36  is 
shifted  or  biased  toward  large  droplet  diameters.  The  data  biasing  is  also  shown  by  the  log- 
probability  plot  of  the  histogram  data  (Fig.  37).  Figure  37  reveals  that  a  simple  log-normal 
distribution  (a  straight  line)  does  not  fit  the  data.  As  in  the  case  of  the  AEDC-PSI  data 
(Section  3.2.2),  the  log-probability  profile  looks  like  a  log-normal  distribution  from  which 
droplets  have  been  removed  (Ref.  17). 

Throughout  the  entire  experimental  program,  biasing  showed  up  in  the  data  obtained 
with  the  commercial-PSI.  Changes  in  test  conditions  did  not  alter  the  bias;  the  output  of  the 
PSI  was  insensitive  to  changes  in  the  spray  nozzle  mass  ratio.  This  is  illustrated  in  Fig.  38,  a 
comparison  between  the  geometric  mean  droplet  diameters  obtained  from  the  PSI  data  and 
these  derived  with  Eq.  (5).  From  the  results  of  Fig.  38,  and  the  experience  gained  with  the 
AEDC-PSI  (Section  3.2.2),  it  was  deduced  that  the  data  biasing  occurred  because  the 
number  density  of  droplets  in  the  simulated  icing  cloud  was  too  large  for  satisfactory 
operation  of  the  PSI. 

A  90-deg  mask  was  added  to  the  collection  lens  of  the  PSI  (Fig.  13)  to  reduce  the 
background  noise  resulting  from  the  large  number  of  droplets  that  cross  the  lasers  outside 
the  crossover  region.  Also,  the  pinhole  in  the  photomultiplier  (PM)  tube  (Fig.  12)  was  made 
smaller  to  reduce  the  effective  size  of  the  sample  volume  viewed  by  the  PSI.  Neither 
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changing  the  pinhole  size  nor  adding  the  mask  reduced  the  data  bias  of  the  PSI  data. 
Consequently,  it  was  concluded  that  the  commercially  developed  PSI  cannot  be  used  in  the 
present  configuration  to  size  droplets  of  a  simulated  icing  cloud. 

3.2.4  BPS  Results 

Figure  39  is  a  typical  histogram  of  the  droplet  size  data  taken  with  the  BPS.  A 
comparison  between  the  data  of  Fig.  39  and  the  holocamera  results  of  Fig.  19  shows  that  the 
BPS  data  are  biased  toward  the  smaller  droplets.  In  fact,  it  is  strikingly  apparent  that  the 
droplet  size  distribution  displayed  in  Fig.  39  is  monodisperse  (i.e.,  95  percent  of  the  droplets 
sized  are  contained  in  a  single  bin  —  the  bin  corresponding  to  a  droplet  diameter  of  8.5  /mi). 

Since  the  accuracy  of  the  holocamera  is  questionable  for  droplet  diameters  smaller  than 
10 /mi,  it  is  reasonable  to  ask  whether  the  distribution  of  Fig.  39  is  not  more  representative  of 
the  conditions  existing  in  the  simulated  icing  cloud  than  is  that  of  Fig.  19.  If  the  droplet  size 
distribution  of  the  cloud  is  monodisperse,  as  shown  in  Fig.  39,  then  the  data  rate  of  the  BPS 
would  be  approximately  500  droplets/sec  (see  Appendix  Q.  However,  during  data 
acquisition  the  average  data  rate  obtained  with  the  BPS  was  approximately  10  droplets/sec. 

Although  a  discrepancy  in  data  rate  hardly  answers  the  question  of  correctness  of  the 
histograms,  it  does  point  out  that  the  BPS  was  not  functioning  as  expected.  A  monodisperse 
panicle  generator  (Ref.  18)  was  used  to  check  the  calibration  of  the  BPS.  This  was  done  by 
producing  a  stream  of  monodisperse  water  droplets  of  known  size  (Ref.  18),  using  the  BPS 
to  obtain  a  measure  of  the  droplet  diameters  of  the  stream,  and  comparing  the  actual 
diameter  sizes  to  the  measured  sizes.  The  results  are  shown  in  Fig.  40,  in  which  the  “line  of 
perfect  agreement”  is  a  step  function  because  the  BPS  data  are  classified  into  bins.  Hence,  a 
droplet  whose  actual  size  falls  within  a  given  bin  width  should  be  sized  by  the  BPS  as  the 
diameter  corresponding  to  the  line  of  perfect  agreement.  Figure  40  shows  that  the  BPS  did 
not  size  the  14  /im  and  20  nm  droplets  correctly.  Each  droplet  was  sized  as  having  an  8.5-/im 
diameter.  These  results  showed  that  the  BPS  indeed  inaccurately  sized  the  droplet  diameters. 
It  was  therefore  concluded  that  the  BPS  could  not  be  used  as  a  droplet  sizer  during  icing 
tests. 


4.0  SUMMARY  AND  CONCLUDING  REMARKS 

A  series  of  experiments  was  conducted  in  the  icing  research  test  cell  to  investigate  the 
particle-sizing  capabilities  of  selected  particle  diagnostic  systems.  The  selected  systems 
include  a  fiber-optics  system,  two  types  of  particle-sizing  interferometers,  and  a 
backscattering  particle-sizing  system.  A  holocamera  was  the  baseline  device  that  provided 
particle  size  data  for  comparison. 
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From  calibration  of  the  holographic  system  it  was  concluded  that  the  smallest  particle 
resolvable  with  the  present  holocamera  is  10  /im.  It  was  also  determined  that  the  calibration 
functions  for  the  holographic  system  change  with  time.  Therefore,  to  ensure  that  the  correct 
calibration  functions  are  used,  repeated  calibrations  are  required  during  data  reduction. 

Calibration  of  the  icing  research  cell  with  the  holocamera  revealed  that  the  diameters  of 
the  droplets  varied  from  10  to  90  /urn  and  that  the  distributions  of  the  droplet  diameters  are 
expressible  as  logarithmic  normal  functions.  Utilizing  the  log-normal  characteristics  of  the 
distributions,  a  central  tendency  and  a  measure  of  the  dispersion  of  the  droplet  diameters 
about  the  central  tendency  were  found  for  each  calibration  condition.  The  mass  median 
diameter  was  the  measure  of  central  tendency  used,  and  the  geometric-standard  deviation 
was  used  as  the  measure  of  dispersion. 

A  correlation  between  the  mass  median  diameter  and  the  spray  nozzle  mass  ratio  is 
presented,  and  a  power  law  was  dervived  to  describe  the  correlation.  The  mass  median 
diameter  decreases  with  an  increase  in  spray  nozzle  mass  ratio;  this  dependence  can  be  used 
as  a  spray  nozzle  calibration.  Thus,  in  the  applicable  range  of  the  calibration  curve,  which 
for  the  mass  median  diameter  is  from  15  to  25  fim,  only  one  droplet  diameter  need  be 
reported  for  each  calibration  condition. 

Data  obtained  with  the  fiber-optics  system  were  also  log-normally  distributed.  The 
power  law  used  that  described  the  calibration  data  also  correlated  the  mass  median 
diameters  obtained  from  the  Fiber-optics  data  and  the  spray  nozzle  mass  ratio.  Although  a 
small  negative  bias  exists  when  the  mass  median  diameters  are  determined  with  the  fiber- 
optics  system,  the  bias  is  within  the  error  limits  of  the  holocamera.  It  is  concluded, 
therefore,  that  the  fiber -optics  systems  can  be  used  in  a  simulated  icing  cloud  environment  to 
determine  the  droplet  size  distribution  and  the  effective  mean  droplet  diameter. 

Since  the  data  obtained  with  both  types  of  particle-sizing  interferometers  were  biased 
toward  large  droplet  diameters,  it  was  concluded  that  neither  the  AEDC-developed  PSI  nor 
the  commercially  developed  PSI  can  be  used,  as  presently  designed,  to  size  droplets  in  a 
simulated  icing  cloud.  The  accuracy  of  the  results  is  dependent  on  factors  related  to  the  high 
number  density  of  water  droplets  rather  than  on  the  inherent  accuracies  of  the  instruments. 
It  was  inferred  that  the  presence  of  large  numbers  of  droplets  within  the  measurement  plane 
increases  the  probability  of  multiparticle  sampling  and  also  severely  increases  the 
background  noise.  These  effects  can  be  minimized  or  eliminated  with  the  use  of  off-axis 
scatter  detection  (Ref.  19).  As  pointed  out  by  Bachalo  et  at.  (Ref.  19),  the  use  of  off-axis 
scatter  detection  extends,  by  an  order  of  magnitude,  the  range  of  particle  number  density  in 
which  an  interferometer  can  operate.  With  such  an  increase  in  the  range  of  operation,  the 
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interferometer  could  become  a  viable  means  of  sizing  droplets  in  a  simulated  icing  cloud.  It 
is  suggested  that  the  off-axis  interferometer  be  evaluated  to  determine  its  effectiveness  in  a 
simulated  icing  cloud. 

The  data  obtained  with  the  backscattering  particle-sizing  system  were  monodisperse  and 
inaccurate.  Thus,  the  backscattering  particle-sizing  system  could  not  be  used  to  measure 
droplet  size  during  icing  tests. 
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Altitude:  4,000  to  22,000  ft 
Horizontal  Extent:  3  miles 
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Pressure  Altitude,  x  10‘ ,  ft 

b.  Intermittent  maximum  (cumuliform  clouds) 
Figure  1.  Concluded. 
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a.  Direct  connect 

Figure  2.  Propulsion  engine  altitude  icing  test  cell  configuration. 
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Figure  2.  Concluded. 
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a.  Schematic  of  general  layout 
Figure  3.  Icing  research  test  cell. 


AEDC-TR-S0-33 


ro 


Viewing  Ports 


Instrumentation  Port 


Connecting 

Ducts 


Plenum  Chamber 


b.  Icing  research  test  cell  looking  upstream 
Figure  3.  Continued. 
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Schematic  icing  spray  system  details 
Figure  3.  Continued. 
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Figure  4.  Research  icing  facility  holocamera. 
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Figure  6.  AEDC  semi-automated  hologram  reconstruction  system. 
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Figure  7.  Fiber-optics  particle-sizing  system. 
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Figure  8.  FOS  optical  arrangement. 
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Figure  9.  FOS  data  handling  system. 
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Figure  11.  AEDC-developed  PSI  optical  system. 
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Figure  12.  Electronics  for  AE DC-developed  RSI 
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Figure  13.  Commercial-PSI  optical  system. 
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Figure  17.  Optical  arrangement  of  backscattering 
particle-sizing  system. 
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a.  Pretest  calibration  of  holographic  system. 


c.  Posttest  calibration  of  holographic  system 
Figure  18.  Holographic  system  calibration. 
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Figure  19.  Typical  histogram  for  research  cell  calibration 
(data  point  No.  3). 
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Figure  21.  Typical  histogram  for  research  cell  calibration 
(data  point  No.  6). 
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Figure  22.  Log-probability  plot  of  holographic  data  point  No.  3. 
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Figure  23.  Log-probability  plot  of  holographic  data  point  No.  9. 
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Figure  24.  Log-probability  plot  of  holographic  data  point  No.  6. 
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Figure  26.  Calculated  LWC  vs  input  LWC. 
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a.  Hologram  number  4,  5.34  cm  from  film  plane 


b.  Hologram  Number  18,  5.66  cm  from  film  plane 
Figure  27.  Reconstructed  holograms. 
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Figure  28.  Typical  FOS  histogram  (data  point  No.  3). 
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Figure  29.  Log-probability  plot  of  FOS  data  point  No.  3. 
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Figure  33.  Typical  log-probability  plot  for  AEDC-PSI. 
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Figure  35.  Comparison  of  AEDC-PSI  data  with  spray  nozzle  calibration. 
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Figure  36.  Typical  histogram  for  commercial-PSI. 
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Figure  37.  Typical  log-probability 


70  80  90  95  98  99  99.9  99.99 

cated  Diameter 


for  commercial-PSI . 


Spray  Nozzle  Mass  Ratio , 

w 

w 

Figure  38.  Comparison  of  commercial-PSI  data  with 
spray  nozzle  calibration. 
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Figure  39.  Typical  histogram  for  BPS. 
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Figure  40.  Measured  vs  actual  droplet  diameters  for  BPS. 
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Table  1.  Icing  Research  Cell  Performance  Capability 


Conditions 

Level 

Maximum 

Minimum 

Pressure  Altitude 

80,000  ft 

Sea  Level 

Inlet  Air  Pressure 

40  psi 

— 

Inlet  Air  Temperature 

185°F 

.  -20°F 

Cell  Mach  Number 

0.7 

— 

Airflow  Rate 

30  lbm/sec 

— 

Atomizing  Air  Pressure 

175  psia 

Atmospheric 

Atomizing  Air  Temperature 

200  °F 

70°F 

Demineralized  Water  Pressure 

290  psia 

Atmospheric 

Demineralized  Water  Temperature 

200°F 

Atmospheric 

Demineralized  Water  Flow  Rate 

12  gal/hr 

0  gal/hr 

Table  2.  Test  Conditions  for  Experimental  Program 


Data 

Point 

Number 

(1/16/78) 

Inlet  Air  Conditions 

Spray  Nozzle  Conditions 

Total 

Temp, 

°F 

Free-Stream 
Bulk  Vel, 
fps 

Alt, 

ft 

Air 

Water 

Press . , 
psia 

Temp, 

°F 

Flow, 

cfm 

Press .  , 
psia 

Flow 

gal/min 

1 

37 

330 

14,500 

112.46 

57 

2.70 

91.67 

67 

0.0462 

2 

38 

330 

14,500 

100.66 

60 

2.30 

89.80 

67 

0.1083 

3 

38 

340 

14,600 

99.72 

58 

2.06 

81  .80 

67 

0.0440 

4 

38 

343 

14,600 

49.90 

67 

1.55 

46.03 

67 

0.0652 

5 

38 

350 

14,600 

25.49 

69 

0.75 

41.64 

67 

0.1874 

6 

39 

335 

14,600 

99.33 

65 

0.93 

105.31 

67 

0.1872 

7 

38 

335 

14,600 

100.03 

67 

2.31 

89.75 

67 

0.1070 

8 

38 

341 

14.600 

125.06 

60 

2.31 

125.89 

67 

0.1861 

9 

38 

339 

14,600 

118-58 

58 

2.55 

103.30 

67 

0.1068 

10 

38 

340 

14,500 

28.62 

68 

0.56 

63.98 

80 

0.1897 

1 1 

38 

343 

14,500 

28.64 

68 

1  .00 

50.63 

80 

0.1175 

12 

38 

339 

14,500 

48,21 

68 

1.82 

89.33 

72 

0.1121 

13 

38 

340 

14,500 

20.36 

68 

1  .07 

30.95 

80 

0.0977 

14 

38 

344 

14,600 

20.33 

68 

0.72 

28.18 

78 

0.0599 
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Table  3.  Parameter  Measurement  Uncertainty 


Parameter 

Range 

Bias  (%) , 

B 

Precision 
Index  (%), 

S 

Degrees 

of 

Freedom 

Uncertainty  (%) , 
±(B  +  t95  S) 

Airflow 

<30  lbm/sec 

0.12 

0.17 

>30 

0.45 

Total  Temperature 

460°R 

0-4 

0.20 

>30 

1,40 

Total  Pressure 

5  to  30  psia 

0.4 

0.30 

>30 

1  .0 

Water  Flow 

1-11.5  gal /hr 

0.4 

0.30 

>30 

1  .0 
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Table  4.  Holocamera  Results 


Data  Point 
Number 

Analytical  Evaluation 

Graphical  Evaluation 

Eq.  (2), 
M, 

Mm 

Eq.  (3), 

S 

g 

Eq.  (4), 
dgm, 
yui 

M, 

ym 

S 

g 

1 

12.6 

1.325 

16.0 

12.5 

1.360 

2 

14.9 

1.387 

20.5 

15.0 

1 .410 

3 

13.3 

1.354 

17.5 

13.2 

1.325 

4 

18.1 

1  .530 

31.3 

18.0 

1 .525 

5 

15.5 

1.471 

24.2 

15-7 

1.490 

6 

21.1 

1.513 

35.4 

21 .0 

1.520 

7 

14.2 

1  .439 

21.1 

14.0 

1.421 

* 

8 

14.3 

1.531 

24.6 

14.2 

1-501 

9 

14.5 

1.403 

20.5 

14.5 

1 .462 

10 

27.5 

1  .465 

42.6 

27.3 

1.400 

11 

24.3 

1  .381 

33.2 

24.0 

1 .362 

12 

18.1 

1  .448 

27.3 

18.0 

1 .450 

13 

12.9 

1  .801 

36.4 

13.0 

1.855 

14 

12.7 

1  .761 

33.1 

12.5 

1.801 
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APPENDIX  A 

METHODS  OF  CALCULATION  FOR  EXPERIMENTAL  DATA 

The  general  methods  and  equations  used  to  compute  the  parameters  presented  in  this 
report  are  given  below.  Where  applicable,  the  arithmetic  average  of  the  pressure  and 
temperature  was  used. 


A-1.0  AIRFLOW 


A-1.1  VENTURI  AIRFLOW 

The  primary  airflows  were  calculated  from  the  equation 


W  =  0.5318 


cr  c,  c  P  A 

I  t  V  s 


(A-l) 


where  C,  =  f(TT,  PT,  Ps),  Cf  =  f(Tr,  PT),  and  Cv  =  f(Ps,  A) 

A-l. 2  ORIFICE  FLOW  (SECONDARY  FLOW) 


W  .  KVA  yfiTp  lPa-Pd) 

The  discharge  coefficient  K  was  defined  by  the  equation 
K.  =  0.5602  +  0.2721/3  -0.7246/92  +  0.8750/83  -  [o.00743  fog  DpipE  -  0.002231  ]  ^  ~  ° 


(A- 2) 


-0.25 

.35 


The  expansion  factor  Y  was  defined  by  the  equation  for  1  >  Pd/Pu  ^  0,63 


Y  =  0.3245  [  —  1 

4  0.6755  4-  [o.033  -0.055/3 

] 

2.703  -/]  -  — 

V  Pu/ 

V  P. 

Pd/Pu  >  0 

/V 

1  r  1 

- 

/P  \ 

Y  =  0.394J  — 

VPu/ 

)-  0.623  4-  0.102  -0.17/3 

1 

- 

_ 

(A-3) 


(A-4) 


(A-5) 
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where  £  —  Dqrifice/Dpipe- 


A-2.0  ALTITUDE 

Altitudes  were  calculated  from  the  tables  of  the  1962  Standard  Atmosphere,  with 


Altitude  =  f  (Ps,  Type  of  Day) 


(A-6) 


A-3.0  VELOCITY  AND  MACH  NUMBER 


The  flight  velocity  and  Mach  number  were  calculated  by  using  an  isentropic  relationship 
with  Pj,  Ps,  and  Ty  and  with  a  ram  recovery  of  unity. 


V 


a 


(A-7) 


(A-8) 


M  = 


V 


a 


(A-9) 


where  Ts  =  f(TT,  Pt,  Ps)- 

A-4.0  FREEJET  VOLUMETRIC  AIRFLOW  RATE 

The  freejet  volumetric  airflow  rate  was  calculated  from  the  relationship 

0  =  A-V 

where  V  =  f  (Tt,  Ps,  Pt). 

A-5.0  LIQUID  WATER  CONTENT 

The  airstream  liquid  water  content  was  calculated  from  the  relationship 


Q  ■  K 

where  Q  =  (A,  Tt,  Ps,  Py)  and  K  =  6.241  x  10" s. 


(A-10) 


(All) 
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APPENDIX  B 

CURVE-FITTING  AND  STATISTICAL  PARAMETERS 
B-1.0  CURVE  FITTING 

The  method  of  least  squares  is  used  to  express  one  variable  in  terms  of  another  even 
though  the  variables  are  not  necessarily  analytical  functions  of  each  other.  In  reporting  the 
data  of  this  report,  two  types  of  least-squares  linear  regressions  were  used. 

B-1.1  LINEAR  FIT 

For  the  data  typical  of  Fig.  B-l,  a  simple  linear  regression  was  used;  the  constants  of  the 
correlating  equation  are  determined  as 


y 


y  =  m  x  +  b 
(Linear  Correlating 
Equation ) 


Figure  B-t.  Typical  Linear  Regression 


n 

□ 

n  £ 

y  x  v 

X  1  7. 

1  j=l '1 

L  i  - 

i  =  l 

n 

T  (2 

‘O’ 

M  x2  Vi=1 

b  = 


(B-l) 


(B-2) 


The  coefficient  of  determination  is  calculated  to  determine  the  quality  of  the  fit  achieved  by 
the  regression.  It  is  given  by 
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Values  of  r2  close  to  unity  indicate  a  better  fit  than  values  close  to  zero. 

B-1.2  POWER  CURVE  FIT 


To  use  the  least-squares  method  to  fit  the  data  of  Fig.  B-2  to  the  nonlinear  power  equation, 
it  was  necessary  to  reduce  the  problem  to  a  linear  one  by  first  taking  the  logarithms,  getting 

fn  y  =  tn  a  +  k  Un  x  (B-4) 

Equation  (B-4)  is  linear  in  the  parameters.  A'  =  fn  a  and  K  =  k.  Hence,  a  linear  regression 
was  applied  to  determine  the  regression  coefficient  of  the  power  curve  fit.  The  coefficients 
are  determined  from 


k  - 


n 


£  ( tn  x.)  2  (  fn  y  , ) 

(fn  x.)  <l'n  y.)  -  — — - - - — - 

1  1  n 


£  (Fn.)2 
1=1 


n 


(B-5) 
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a  = 


f  * 


exp  < 


L 


(B-6) 


The  coefficient  of  determination  is  given  by 


(B-7) 


B-2.0  STATISTICAL  PARAMETERS 

Let  a  set  of  data  (xj,  yj)  be  represented  by  the  least-squares  curve  fit  equation 

y  =  f(x)  (B-8) 


Also,  let 


®i  =  f(*j>  "  Tj  (B-9) 

be  the  deviations  of  the  computed  values  from  the  observed  values.  Then  the  mean 
deviation,  or  error,  is  defined  as 


S  = 


,?,6' 


(B-10) 


and  the  standard  deviation  of  the  errors  is  given  by 
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APPENDIX  C 

CALCULATED  DATA  RATE  FOR  BACKSCATTERING 
PARTICULE  SIZING  SYSTEM 

The  ayerage  liquid  water  content  of  an  icing  cloud  is  defined  as 

LWC  =  -  pd3  .  N  (C-l) 

6 


where  d  is  some  average  droplet  size  of  the  icing  cloud.  If  the  cloud  has  a  monodisperse 
droplet  distribution  as  reported  for  the  backscattering  particle-sizing  system  described  in 
Section  3.2.4,  then  the  average  LWC  becomes  the  actual  LWC.  Hence,  to  determine  the 
droplet  number  density  for  a  monodisperse  distribution, 

N  =  LWC 


For  the  histogram  of  Fig.  39,  d  =  8.5  pm,  input  LWC  =  0.35  g/m3;  hence, 


) 


(C-2) 


N  =  0.35 


1,088  particles/ cc 


The  probe  volume,  Vs,  is  0.5  mm3  and  the  firing  frequency  of  the  laser,  q  is  lOVsec  (Section 
2.2.4).  The  data  rate  of  the  backscattering  particle-sizing  system  becomes 

N  =  V>N  (C-3) 


or  for  the  histogram  of  Fig,  39 


N  =  (0.5  x  L0‘3)  (103)  (1,088) 


=  504.4  particles/ sec 


> 
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NOMENCLATURE 

A  Area,  ft2,  in.2 

a  Constant  or  sonic  velocity,  ft/sec 

B  Bin  number  for  FOS 

b  Constant 

Cd,  Cv  Flow  coefficients 

Cp  Specific  heat  at  constant  pressure,  ft2/sec2°R 
Q  Thermal  expansion  coefficient,  in. /in. 

D0  Orifice  diameter,  ft,  in. 

Dp  Pipe  diameter,  ft,  in. 

d  Mean  effective  droplet  diameter,  /mi 

da  Actual  diameter,  /tm 

dgm  Mass  geometric  mean  diameter,  /im 

d,  ith  droplet  diameter,  /*m 

dm  Measured  diameter,  /im 

fi  Frequency  of  observation  of  diameter  dj 

FOS  Fiber-Optics  System 

gc  Dimensional  constant,  lbm-ft/lbf-sec2 

K  Discharge  coefficient 

k  Constant  or  conversion  factor 

LWC  Liquid  water  content,  gm/m3 

M  Geometric  mean  droplet  size  based  on  count,  /im,  or  Mach  number 
m  Constant 

N  Droplet  number  density,  particles/cc 
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N  Data  rate,  samples/sec 

n  Sum  of  all  droplets  or  sum  of  all  data  points 

nj  Number  of  droplets  contained  in  the  ith  line 

Pa  Atomizing  air  pressure,  psfa 

Pd  Static  pressure  downstream  of  orifice,  psfa 

Ps  Static  pressure,  psfa 

Pt  Total  pressure,  psfa 

t 

Pu  Static  pressure  upstream  of  orifice,  psfa 

Pw  Spray  nozzle  water  pressure,  psfa 

Q  Volumetric  flow  rate,  ft3 /sec 

r2  Coefficient  of  determination 

Sg  Geometric  standard  deviation 

Sj  Standard  deviation  of  measurement  errors 

Tt  Total  temperature,  °R 

V  Air  velocity,  ft/sec 

W  Primary-airflow  rate,  Ibm/sec 

Wh2o  Spray  nozzle  water  flow  rate,  lbm/sec 
wa  Mass  flow  rate  of  spray  nozzle  atomizing  air,  lbm/sec 

w»  Mass  flow  rate  of  spray  nozzle  water,  lbm/sec 

Y  Expansion  factor 

0  D0/Dp 

7  Ratio  of  specific  heats 

ith  deviation  between  the  computed  and  observed  values 
e  Density,  lbm/ft3 
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